Thermal broadening of the J-band in disordered linear molecular aggregates: 

A theoretical study 
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We theoretically study the temperature dependence of the J-band width in disordered linear 
molecular aggregates, caused by dephasing of the exciton states due to scattering on vibrations of the 
host matrix. In particular, we consider inelastic one- and two-phonon scattering between different 
exciton states (energy-relaxation-induced dephasing), as well as elastic two-phonon scattering of the 
excitons (pure dephasing). The exciton states follow from numerical diagonalization of a Frenkel 
Hamiltonian with diagonal disorder; the scattering rates between them are obtained using the Fermi 
Golden Rule. A Debye-like model for the one- and two-phonon spectral densities is used in the 
calculations. We find that, owing to the disorder, the dephasing rates of the individual exciton 
states are distributed over a wide range of values. We also demonstrate that the dominant channel 
of two-phonon scattering is not the elastic one, as is often tacitly assumed, but rather comes from 
a similar two-phonon inelastic scattering process. In order to study the temperature dependence 
of the J-band width, we simulate the absorption spectrum, accounting for the dephasing induced 
broadening of the exciton states. We find a power-law (T p ) temperature scaling of the effective 
homogeneous width, with an exponent p that depends on the shape of the spectral density of host 
vibrations. In particular, for a Debye model of vibrations, we find p ~ 4, which is in good agreement 
with experimental data on J-aggregates of pseudoisocyanine [J. Phys. Chem. A 101, 7977 (1997)]. 
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PACS numbers: 71.35.Aa; 36.20.Kd; 78.30.Ly 



I. INTRODUCTION 



Ever since the discovery of aggregation of cyanine dye 
molecules by Jelleyi and Scheibe^ the width of the ab- 
sorption band of these linear aggregates (the J-band) has 
attracted much attention. At low temperatures, the J- 
band may be as narrow as a few tens of cm 's (for pseu- 
doisocyanine), while at room temperature it typically is 
a few hundred cm" 1 ^ The small width at low temper- 
ature is generally understood as resulting from the ex- 
citonic nature of the optical excitations^ which leads 
to exchange narrowing of the inhomogeneous broaden- 
ing of the transitions of individual molecules^ As the 
optically dominant exciton states in inhomogeneous J- 
aggregates occur below the exciton band edge, their 
vibration-induced dephasing is strongly suppressed at 
low temperatures. Thus, the J-band is inhomogeneously 
broadened, except for a small residual homogeneous com- 
ponent (~ 0.1 cm -1 ) caused by spontaneous emission of 
the individual exciton states underlying the spectrum. 
Upon increasing the temperature, the vibration-induced 
dephasing of the exciton states increases, the J-band 
broadens, and obtains a more homogeneous character. 

Over the past twenty years, the temperature depen- 
dence of the J-band width and the homogeneous broad- 
ening of the exciton states, has been studied by several 
authors. In 1987, De Boer, Vink, and WiersmaS, per- 
formed accumulated-echo experiments to study the tem- 
perature dependence of the pure dephasing time in J- 
aggregates of pseudo-isocyanine bromide (PIC-Br) and 
found that within the temperature range 1.5 K to 100 
K, this time is linearly proportional to the occupation 



number of a mode with an energy of 9 cm . They as- 
signed this to a librational mode of the aggregate. Later 
on, Fidder, Knoester, and Wiersma- (see also Ref. 
and 0) showed that a similar study carried our over a 
wider temperature range, 1.5 K to 190 K, required the 
occupation numbers of three vibration modes, at 9 cm , 
305 cm -1 , and 973 cm" 1 . Using the hole burning tech- 
nique, Hirschmann and Friedrichp* studied the homoge- 
neous width of the exciton states in pseudoisocyanine 
iodide (PIC-I) over the temperature range 350 mK to 80 
K. They were able to fit their measurements by a super- 
position of two exponentials, with activation energies 27 
cm" 1 and 330 cm" 1 , and they attributed the broadening 
to scattering of the excitons on an acoustic mode and an 
optical mode, respectively, of the aggregate. Finally, in 
1997 Renge and Wildii measured the temperature de- 
pendent width, A(T), of the total J-band of pseudoiso- 
cyanine chloride (PIC-C1) and fluoride (PIC-F) over the 
wide temperature range 10 K to 300 K. They found that 
over this entire range their data closely obeyed a power- 
law dependence A(T) = A(0) 4- bT p with the exponent 
p = 3.4. They suggested the scattering of the excitons 
on host vibrations as a possible source of this behavior. 

The above overview clearly demonstrates that the tem- 
perature dependence of the most important characteristic 
of J-aggregates, namely the J-band, is not understood. It 
even is not clear what is the source of thermal broaden- 
ing: dephasing due to vibrational modes of the aggre- 
gate itself or due to modes of the host matrix. Theoreti- 
cally, the study of the exciton dephasing in J-aggregates 
is complicated by the fact that static disorder plays an 
important role in these systems, as is clear from the 
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strongly asymmetric low-temperature lineshapeA* The 
simultaneous treatment of disorder, leading to exciton 
localization, and scattering on vibrational modes, is a 
problem that requires extensive numerical simulations. 
Such simulations have been used previously to model the 
optical response of polysilanesjiS*^ the temperature de- 
pendent fluorescence 14 and transport 1 ^ properties of J- 
aggregates, and single molecular spectroscopy of circular 
aggregates^ Alternatively, stochastic models have been 
applied to describe the scattering of excitons on vibra- 
tions in disordered aggregates and the resulting spectral 
line shap o 17 i 18 and relaxation dynamics^ 

In this paper, we report on a systematic theoretical 
study of dephasing of weakly localized Frenkel excitons 
in one-dimensional systems, focusing on the effect of 
the scattering of the excitons on the vibrational modes 
(phonons) of the host. For a chain-like configuration it 
seems physically reasonable to assume that the coupling 
to host vibrations dominates the dynamics of the exci- 
tons. For coupling to vibrations in the chain, one ex- 
pects self-trapped exciton states^ for which in most ag- 
gregates no clear signature is found. We describe the 
phonons by a Debye model and consider one-phonon as 
well as (elastic and inelastic) two-phonon contributions to 
the dephasing. The scattering rates between the various, 
numerically obtained, exciton states are derived using the 
Fermi Golden Rule. Due to the disorder, the dephasing 
rates of individual exciton states are distributed over a 
wide range, in particular at low temperature, making it 
meaningless to associate the width of individual exciton 
levels with the J-band width. Rather, this width is deter- 
mined by direct simulation of the total absorption spec- 
trum and is found to scale with temperature according 
to a power law. 

The outline of this paper is as follows: In Sec. [H] 
we introduce the Hamiltonian for excitons in a disor- 
dered chain and coupled to host vibrations. General ex- 
pressions for the exciton dephasing rates are derived in 
Sec. IIIII Section IIVI deals with the temperature depen- 
dence of the dephasing rates for the disorder-free case, 
where analytical expressions can be obtained. In Sec. 
we give the results of our numerical simulations for the 
dephasing rates in the presence of disorder and analyze 
their temperature dependence and fluctuations as well 
as the total J-band width. We compare to experiment in 
Sec. lVII and discuss an alternative mechanism of dephas- 
ing due to the coupling of excitons to a local vibration. 
In Section IVIII we present our conclusions. 



II. MODEL 

We consider an ensemble of J-aggregates embedded in 
a disordered host matrix. The aggregates are assumed to 
be decoupled from each other, while they interact with 
the host. A single aggregate is modeled as an open lin- 
ear chain of N coupled two-level monomers with parallel 
transition dipoles. The interaction between a particular 



monomer and the surrounding host molecules in the equi- 
librium configuration leads to shifts in the monomer's 
transition energy. Due to the host's structural disor- 
der, this shift is different for each monomer in the ag- 
gregate, giving rise to on-site (diagonal) disorder. More- 
over, vibrations of the host couple to the aggregate ex- 
cited states, because the associated displacements away 
from the equilibrium configuration dynamically affect the 
monomer transition energies. Accounting for these shifts 
up to second order in the molecular displacements, the 
resulting Hamiltonian in the site representation reads: 

H = H C * + H hath + V {1) + l/ (2) , (la) 

with 

N N 
HC * = Y. £ n\n)(n\+ Jnm\n)(m\, (lb) 

n— 1 n,rn—l 

^bath = £ ^ a t a? ( (lc) 
9 

JV 

1/(1) = EE f 4 1) whk+4) ! ( ld ) 

n— 1 q 

N 

y(2U EE 1/ Si ?i >HK+ fl l)K+4)- ( le ) 

n— 1 qq' 

Here, £P X is the bare Frenkel exciton Hamiltonian, with 
|n) denoting the state in which the nth monomer is ex- 
cited and all the other monomers are in the ground state. 
The monomer excitation energies, £\, e%, . . . , ejy, are un- 
corrected stochastic gaussian variables, with mean e and 
standard deviation a, referred to as the disorder strength. 
Hereafter, e is set to zero. The resonant interactions J nm 
are considered to be nonrandom and are assumed to be 
of dipolar origin: J mn — —J/\n — m\ 3 (J nn = 0), with 
J > denoting the nearest-neighbor coupling. 

i? bath describes the vibrational modes of the host, la- 
beled q and with the spectrum ui q (h = 1). The oper- 
ator a q (a q ) annihilates (creates) a vibrational quantum 

in mode q. Finally, the operators V r< - 1 - ) and describe 
the linear and quadratic exciton-vibration coupling, re- 
spectively, where the quantities Vnq and V^ q , indicate 
their strengths. We do not derive explicit expressions for 
these coupling strengths, as we aim to treat them on a 
phenomenological basis. In particular, owing to the dis- 
ordered nature of the host, we consider these strengths 
stochastic quantities, for which we only specify the fol- 
lowing stochastic properties with respect to the site index 
n: 
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v mq v nq > — °n 



rnqq' nqq' 



yd) 



V {2 ) 

QQ 



(2b) 



(2c) 



where the angular brackets denote averaging over real- 
izations of V$ and V^ ql . The Kronecker symbol in 
Eqs. I)2bjl and J2cJ implies that the surroundings of dif- 
ferent monomers in the aggregate are not correlated. 

In general, the Hamiltonian Eq. can not be diago- 
nalized analytically. If the exciton-vibration coupling is 
not too strong, the method of choice is first to find the ex- 
citon eigenstates by numerical diagonalization and then 
to consider the scattering on the basis of these eigen- 
states. Explicitly, the exciton states follow from the 
eigenvalue problem 



JV 

E 



rrex 



i/=l,2... 



(3a) 



helium temperature and becomes several times broader 
at room temperature, and (ii) - the absence of a fluo- 
rescence Stokes shift of the J-band (see, e.g., Ref. 0). 
The extended nature of the exciton states in J-aggregates 
helps to reduce the exciton-vibration coupling, as it leads 
to averaging of the static as well as dynamic fluctua- 
tions of the site energies, effects known as exchange^ and 
motional2iS narrowing, respectively. The weakness of 
the exciton-vibration coupling allows one to calculate the 
scattering and dephasing rates of the excitons through 
perturbation theory. This analysis is presented in the 
next section. 



III. DEPHASING RATES 

Following the arguments given at the end of Sec. [HI 
we will use Fermi's Golden Rule to calculate the rate for 
scattering of excitons from one localized state, \v), to 
another one, The result reads 



(nJ\v®\v,i) 



where HI 



(n\H c 



and E v is the eigenenergy of 



the exciton state 



N 



n=l 



(3b) 



In the exciton representation, the exciton-vibration in- 
teractions take the form 



JV 



v(1) = E ( 4s 

H,U=l q 



x 5(E M 



a, 



(6) 



Here, the superscript £ = 1,2 distinguishes between one- 
and two-vibration-assisted exciton scattering. Further- 
more, fli and £lf are the energies of the vibration bath in 
the initial = \{n q }i)) and final (|/) = \{n q }f)) states, 
respectively, where {n q } denotes the set of occupation 
numbers of the vibrational modes. The quantity p(f2i) is 
the equilibrium density matrix of the initial state of the 
bath. Finally, the angular brackets indicate that we av- 
erage over the stochastic realizations of the surroundings 
of each monomer in the aggregate. 



N 

y{2) = E E^Im)(hk+4)K'+4') ! ( 4b ) 

fi,u—l qq' 

where Vffiq and V^J qq , are the matrix elements of the 
vibration- induced scattering of an exciton from state \v) 
to state given by 



N 



v {1) - V 

V y 9 _ / i nq Vim 



T ,(2) _ T/ (2) 
'livqq' ~ /-^l V nqq' l P^n i Pun ■ 
n=l 



y 



(5a) 



(5b) 



Spectroscopic data on J-aggregates clearly reveal that 
the exciton-vibration coupling in these systems is usu- 
ally weak. For the prototypical J-aggregates of PIC, this 
claim is corroborated by two facts: (i) - the narrowness 
of the J-band, which only is a few tens of cm -1 at liquid 



A. Linear exciton-vibration coupling 

In a one-phonon- assisted scattering process, the oc- 
cupation number of one phonon mode q increases or de- 
creases by one, corresponding to emission and absorption 
of a vibrational quantum, respectively. Consequently, 
Substituting the explicit form of the 



±Ulq. 



operator V^ 1 ) from Eq. I|ld|l into Eq. (JBJ, we obtain 



JV 



W« = 2^5»LEK (1 

n—1 q 

x [[n(u; g ) + 1] Siuifiv +w 9 ) 
+ n{uj q ) 5(w fiU - uj q ) 



(7) 



where n(u> q ) = [exp(oj g /T) — is the mean occupation 
number of the vibrational mode q (the Boltzmann con- 
stant fcjg = 1) and iti^y = — E v . In deriving Eq. J7J, 
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we used the properties of the stochastic function v£q 
given by Eqs. (|2a|) and l|2b(l . Defining the one- vibration 
spectral density as 



F^(cj) = 2TrJ2\v q {1) \ 5{u-u q ) 
i 

we can rewrite Eq. (0 in the form 

N 



(8) 



n=l 



n(—u)p V ) + 1, uj^ < 



(9) 



As we observe, the rate W$ is proportional to the 
overlap integral of the site occupation probabilities, tp^ n 
and ifi1 n , of the exciton states involved. First of all, this 
leads to a strong suppression of the scattering rate if 
states \n) and \v) overlap weakly or not at all. Second, 
as the low-energy exciton states in a disordered chain 
exhibit large fluctuations in their localization size^ also 
the scattering rates may undergo large fluctuations (see 
Sec.©. 

The dependence of W$) on the energy mismatch w M „ is 
determined by the one-phonon spectral density J-^iuj). 
Characterizing this function requires knowledge of the 
vibrational spectrum u> q as well as the q dependence of 
the exciton- vibration coupling Vq X \ For the special case 
of scattering on acoustic phonons of a relatively long 
wavelength, we have ^^(lu) ~ lu 3 . This behavior re- 
sults from the ^-dependence of the density of states of 
acoustic phonons, combined with the fact that in the 



long- wavelength limit \V q 



(1)|2 



^,2^25 Q ne ma y con _ 



sider this a Debye-like model, in which one replaces the 
summation over the mode index q by an integration over 
the frequency tu q according to the well-known rule: 



E 



c 



(10) 



Here, C is an irrelevant constant, which we will incorpo- 
rate in an overall free parameter (see below), and lu c is 
a cutoff frequency. It is important to note that uj c is not 
necessarily related to the Debye frequency: the generally 
very complex density of vibrational states in a disordered 
solid may on average exhibit an uj 2 scaling up to a given 
frequency u> c . 

Inspired by the above, we consider a slightly wider class 
of one-phonon spectral densities, given by 

^(u)=W^(jY &{u c -u). (11) 

Here, W^ 1] is a free parameter in the model, which 
characterizes the overall strength of the one-vibration- 
assisted scattering rates. It absorbs a number of other 



parameters characteristic for the host lattice (such as 
the velocity of sound), the constant C from Eq. (TDJ, 
as well as the strength of the transfer interaction J (for 
details, see Ref . I25I) . <d(x) is the Heaviside step function. 
When performing numerical simulations, we will mostly 
use a = 3, for which the spectral density of acoustic 
phonons in the long-wave limit is recovered. In some in- 
stances, however, we will discuss how results depend on 
the exponent a. 

Several observations support considering a Debye- 
like vibration spectral density, even for a disordered 
host. Thus, for strongly disordered Yb 3+ doped 
phosphate glasses, a parabolic behavior of the one- 
phonon spectral density was found over a rather wide 
range of measurement (0 to 100 cm -1 )i2& Furthermore, 
closely related spectral densities of the form !F^(oj) ~ 
(uj/uj c ) a exp(— u/u) c ) (or linear combinations of such 
functions) have been used successfully to fit the optical 
dyna mics in phot osynthetic antenna complexes (see, e.g., 
Refs. l27l28l2<l30D . 

One-phonon-assisted scattering results in the transi- 
tion of an exciton from a given state \u) to state |//), 
where necessarily fx ^ v. In other words, this type of 
scattering changes the occupation probabilities of the ex- 
citon states and thus causes population (or energy) relax- 
ation. The population relaxation in turn contributes to 
the dephasing of state \v). The corresponding dephasing 
rate is given by (see, e.g., Ref. l3l|) : 



r« = - 

v ~ 2 



(12) 



Thus, represents the one-phonon-assisted contribu- 
tion to the homogeneous broadening of the excitonic level 
v. We note that the indirectly depend on tempera- 
ture through the n{LO^ v ) [cf. Eq. @]. The temperature 
dependence of the sum over scattering rates in Eq. H12|) 
and the corresponding width of the total exciton absorp- 
tion spectrum will be analyzed in Sees. HVl andlVl 



B. Quadratic exciton-phonon coupling 



When excitons scatter on the second-order displace- 
ments of the host molecules, described by the operator 
, the occupation numbers of two phonon modes q 
and q' with frequencies to q and u) q i change by ±1. Thus, 
VLf — fij = ±uj q ± w q i, where any combination of plus 
and minus is allowed. The corresponding scattering rates 

("2) m 

W^J are obtained from Eq. taking into account the 
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stochastic properties of V^ q , given by Eqs. (|2a|l and ({2"cj) 

N 2 

= 2t: Yip 2 up 2 Y\V (2 ) 

n—1 qq' 

x [[n(w g ) + l] [n{uo q i) + l] 8(uo^ + uj q + bj q >) 

+ 2n(L0q)[n(L0q') + 1]S(UJ^ -UJq+LOq') 



+ n(LU q )n(uJqi) 5(UJ^ - UJq - UJ q >) 



(13) 



If in analogy to the one-vibration-assisted scattering, we 
define the two- vibration spectral density T 1 ^ (lj, lj') as 

^{w,w') = 2irY d \vS ] \ ^-io q )5{w' , (14) 

qq' 

(2) 

the scattering rate W^J takes the form 

N 

71 = 1 J 

x [ [n(uj) + l] [n(uj') + 1] (5(w M „ + lj + lj') 
+ 2n{uj) [n(u)') + l] S(lj^ - lj + lj') 



+ n(u))n(u>') 5(0)^ — lj — lj') 



(15) 



Similar to T^(uj) [Eq. ©], we will use a parametrization 



Wn (2) fLJLJ 



J \ J 2 



e^c-^e^c-w') , (i6) 



(2) 

where Wj is a free parameter that characterizes the 
overall strength of the two-vibration-assisted scattering 
rates [cf. and the discussion following Eq. lfTT)l], 

Two main types of two-phonon-assisted processes may 
be distinguished. Similarly to the one-phonon case, an 
inelastic channel exists, where scattering occurs between 
different exciton states, thus giving rise to population 
relaxation. However, also an elastic channel is present, in 
which an exciton is scattered by emitting and absorbing 
a phonon of the same energy, and the final exciton state 
is identical to the initial one. This process results in pure 
dephasing of the exciton state, with a rate given by 

N 

W^=2j2ft n / d^(,,,)n(,)[n(,) + l] . (17) 



The quantity Yl n =x ftn 1S recognized as the inverse par- 
ticipation ratio, 3 - 2 which is inversely proportional to the 
localization size of the exciton state \v). Thus, we see 



that pure dephasing is suppressed for more extended 
states, an effect known as the motional narrowingi^iiSS 

Like in the one-phonon assisted process, the scatter- 
ing rate between different states \v) and \/jl) is propor- 
tional to the overlap of their site occupations. The final 
expressions for the rates W/£ resulting from Eqs. H15|) 
and l|16fl are derived in Appendix A; they depend on the 
sign of LJ^ as well as on the relation between {lj^I and 
lj c . Distinction is made between three inelastic channels: 
downward (||), in which two phonons are emitted, cross 
(Ti), m which one phonon is absorbed and another is 
emitted, and upward (ft), in which two phonons are ab- 
sorbed. The fourth type of scattering is the elastic (pure 
dephasing) channel, discussed above already. When cal- 

(2) 

culating the two-phonon assisted dephasing rate Tl ' of 
state we will account for elastic as well as inelastic 
contributions: 



r< 2 > = - w {2) 

v 2 " ' 



(18) 



This rate depends on temperature as a result of the mean 
occupation numbers h(uj) and n(u/) of the vibrational 
modes. 



IV. DISORDER-FREE AGGREGATE 

In order to gain insight in the temperature dependence 
of the dephasing rates and the absorption bandwidth, it is 
useful to start by considering a homogeneous aggregate, 
i.e., a = 0. Analytical results can then be obtained if we 
restrict the resonant interactions J nm to nearest-neighbor 
ones. In this approximation (which we will relax in our 
numerical analysis), we have 



N + l 



1/2 



ixvn 



N+l ' 



E„ 



-2 J cos ■ 



N+l 



(19a) 



(19b) 



The corresponding overlap integrals occurring in Eqs. Q 
and Ijlijfl now read 



N 



E2 2 



1 



N + l 



1+2 OW + &lt+v>N+l) 



(20) 



A. One-phonon-assisted dephasing 

In a homogeneous linear chain, the lowest exciton 
state, \v = 1), contains almost all oscillator strength, 
thus dominating the absorption spectrum»Sa£ Therefore, 
the dephasing rate of this state, , is of primary inter- 
est. As follows from Eq. (|12fl . it is determined by the sum 
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over scattering rates to the other exciton states (/i ^ 1), 
all of which are higher in energy. 

In order to evaluate I^ , we replace the summation in 
Eq. U2J) by an integration, (JT -> [(N + l)/n]JdK), 
which is allowed if TV > 1 and T > E 2 - E\. We will 
also assume that T <C J, which implies that the relevant 
exciton levels are those near the lower exciton band edge, 
where = -2J+JK 2 with K = 7r/x/(jV+l). Changing 
the integration variable to x = JK 2 /T, using E^ — E\ s» 
JK 2 and replacing the lower integration limit ir/(N + 1) 
by zero, we obtain 



.(i) _ w 



(i) 



4-7T 



cix ■ 



1 



(21) 



For temperatures T<Cw Cl the upper integration limit 
may be extended to infinity, and we arrive at 



.(i) 



(i) 



4/T 



c «+ 



1 



a+4 



(22) 



where T(z) and C(z) are the gamma-function and the 
Rieman zeta-function, respectively. Thus, for T « u c 
the one-phonon-assisted dephasing rate shows a power- 
law temperature dependence. Note that for our model 
of acoustic phonons (a = 3), increases quite steeply, 

namely as T 7 / 2 . From numerical evaluation of 1^ for 
a homogeneous chain with all dipole-dipolc interactions, 
we have found that the exponent 7/2 is increased to 3.85, 
mainly as a consequence of logarithmic corrections in the 
exciton dispersion near the lower band edge^*^ If we go 
beyond the parabolic range of the energy spectrum, the 
growth becomes even steeper; the exponent then tends 
to 4, because the density of states becomes a constant 
towards the center of the band. 

In the opposite limit T > w c , the exponential in the 
denominator of Eq. (|21|l can be expanded in a Taylor 
series. Up to second order, one obtains 



(i) 



2?r(2a - 



-f-r 



T 
1 



(23) 



which simply reflects the linear high-temperature depen- 
dence of the mean occupation number n(u>). Obviously, 
this scaling also holds in the presence of disorder. 



B. Pure dephasing 

We now turn to the temperature dependence of the 

(2^ (2^ 

pure dephasing rate T„ = (1/2) Wuv . As for a homoge- 
neous aggregate this rate does not depend on the state 

(2) 

index v, we will simply denote it as r purc . Using the ex- 
plicit form of (u>, lo') given by Eq. (|17fl . we arrive at 



From Eq. QZQ it follows that for T < lo c (lo c /T -> oo), 



p(2) 

pure 



3 W, 



(2) 







2N + 1 



r(2a- 



l.K(2ni(^ 



2a+l 



while for T > u c (w c /T->0), 



r(2) _ 3 m (2) 

purc 2 (2a- 1)(JV + 1) 



jt(t) 



2a-l /J 1 
J 



(25a) 



(25b) 



For the case of scattering on acoustic phonons (a = 3) 

and T « w Cl we thus arrive at Tpure oc T 7 . This temper- 
ature dependence resembles that for the pure dephasing 
of an isolated state of a point center, derived by McCum- 
ber and Sturge^ The only difference is that the exciton 
dephasing rate undergoes suppression by a factor of N+l 
due to the motional narrowing effect. We note that this 
narrowing is not observed for one-phonon-assisted de- 
phasing [see Eqs. JUJ) and ffity]. The T 2 -scaling of 
in the high-temperature limit (T lu c ) results from the 
square of the mean phonon occupation number involved 
in Eq. fT7|. 

To conclude this section, we stress that the T 2a+1 and 
T 2 scaling relations of the pure dephasing rate with tem- 
perature obtained here, also hold for disordered aggre- 
gates, because this result is determined only by the two- 
vibration spectral density !F^ 2 '(u),u>'). The suppression 
factor, however, will then be determined by the exciton 
localization size, rather than the chain length. 



C. Inelastic two-phonon-assisted dephasing 

Finally, we analyze the dephasing rate of the super- 
radiant state (y = 1) resulting from the two-phonon in- 
elastic scattering of excitons. This rate, which will be 

(2) 

denoted as is determined by the sum of scattering 



rates to all higher states, 1^ = (1/2) J2^i W^ 1 . Using 
Eqs. (|15l) and lltj|) . and making the same assumptions as 
in the case of one-phonon-assisted dephasing (Sec. lIV A"j l. 
we obtain 



r(2) 



.(2) 



W (2) ( T 



4tt \J 

x a y a 
(e x - l)(ev - 1) 

2e v 



2a+f r u c /T i-uJc/T 

dx / dy 
(i Jo 



9(x-y) + 



1 



V s + y 



In the low-temperature limit, T <C lo c {uj c /T 
Eq. yields 



(26) 



oo 



3 W {2) 
p(2) = 6 W 

pure 2N+1 



2a+l 



dx 



(e* - If 



(24) 



(2) = kW£ 
incl 4tt \J 



(27) 
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where the numerical factor k is given by the double in- 
tegral in Eq. I|26l) . with both upper limits replaced by 
infinity. Comparing Eq. (|27Jl with Eq. (|25a|l . we see that 
I-inci is characterized by a steeper temperature depen- 

(2) 

dence than Tpurc- In particular, for a = 3 the rate 

■^inei oc T 15 / 2 . This means that at higher tempera- 
tures the inelastic two-phonon channel of dcphasing can 
compete with the elastic one (Sec. IV B|) . In the high- 
temperature limit, T ^> uj c {oj c /T — > 0), the rate 
is proportional to T 2 , which is the same scaling relation 
as in the case of pure dephasing [Eq. Ij25b|l ]. 



V. DISORDERED AGGREGATES 

In this section, we will analyze the temperature de- 
pendence of the various dephasing contributions in the 
presence of disorder and their effect on the width of 
the J-band. As we will see (and already anticipated in 
Sec. IIII A|l , the dephasing rates are spread over a wide 
region, in particular at low temperature. As a conse- 
quence, it is not clear a priori what value (mean, typical, 
or other) of these rates should be related to the homoge- 
neous width of the absorption spectrum. Therefore, the 
width of the J-band was obtained by direct simulation of 
the absorption spectrum, using the calculated dephasing 
rates to broaden each of the exciton transitions in the 
band. Explicitly, we have 



ME) 



(E„=i <p» 



7T (E — E v 



r 2 



(28) 



is the dimensionless oscillator 



where F v 

strength of the uth exciton state and Y v = 7„/2 + T 

(2) 

T;> ' is the total homogeneous width of this state. Here, 
rW and I^ 2 ) are given by Eqs. ifT^ and (jT%|> . respec- 
tively, and 7„ = ^qF^ is the radiative decay rate of state 
v (70 denotes the radiative constant of a monomer). In 
all the simulations, we will assume the limit lu c 3> T, 
which implies that u) c — > 00 in Eqs. and l|16|l . As 
before, the angular brackets denote the average over the 
random realizations of the site energies {£ n }- The re- 
sulting J-bandwidth A was determined as the full width 
at half maximum of the thus calculated absorption spec- 
trum. A similar approach has been used2& to simulate the 
absorption spectrum of THIATS aggregates^ at room 
temperature. 



A. One-phonon-assisted dephasing 

In order to study fluctuations in the dephasing rates, 
we analyzed their statistics, focusing on the lowest exci- 
ton state for each randomly generated disorder realiza- 
tion. This choice was motivated by the fact that the 
low-lying states dominate the absorption spectrum. The 



dephasing rate of the lowest state is denoted rl 1 -*; its 
distribution, collected by considering 3 x 10 4 disorder re- 
alizations for chains of N = 500 molecules and a disorder 
strength a = 0.135 J is presented in Figs.^a) and^Jb). 
In generating these figures, we used a one- vibration spec- 
tral density of the form Eq. (fill) , with a = 3. Further- 
more, the rates were calculated for T — 0.06 J [Fig.^a)] 
and T = 0.23 J [Fig. ^b)]; for the protoptypical aggre- 
gates of pseudoisocyanine (J = 600 cm -1 ), this agrees 
with temperatures of 50 K and 200 K, respectively. We 
note that the horizontal axis of the distributions is scaled 
by the average value r^, so that the value of W^ 1 does 
not affect the figures. Of course, this scaling renders the 

axis temperature dependent, because strongly de- 
pends on T, as we will see below (Fig. |3J). 

These figures clearly demonstrate that the relative 
spread in may be considerable, in particular at 

low temperatures. This may be understood from the 
local band-edge level structure of the disordered tight- 
binding Hamiltonian^ For temperatures smaller than 
the J-bandwidth, the exciton scatters between discrete 
levels in the vicinity of the band edge that are local- 
ized in the same region of the chain. Both the en- 
ergy spacing between these states and their localiza- 
tion size undergo large fluctuations: SE^ ~ and 



S (T,n=i vlnVln) ~ En=i flnfln^ As a consequence, 



sw. 



(1) 



- WjhJ . When the temperature is increased, the 
exciton in the lowest state may scatter to many higher- 
lying states, which often are delocalized over an apprecia- 
ble part of the chain. This smears the fluctuations that 
occur in the scattering rates between individual states, 
leading to a decrease in the relative spread of the de- 
phasing rate of the lowest state. 
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F (2) /r (2) 
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FIG. 1: Normalized distributions of the one-phonon-assisted 
and two-phonon-assisted dephasing rates ri [(a) and (b)], 
[(c) and (d)], and T?) [(e) and (f)] of the lowest exciton 
state in each one of 3 x 10 4 disorder realizations for chains of 
500 molecules. Upper panels correspond to T = 0.06J, while 
lower panels were calculated for T = 0.23J. All plots were 
obtained for a disorder strength a — 0.135 J and a one- or 
two-vibration spectral density given by Eq. JTTJ or Eq. JTgJ 
with a = 3. Note that the horizontal axes are scaled by the 

average values of the three rates considered, indicated as , 
(2) (2) 

rff , and r Tl \ These averages depend on temperature. 
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0.1 0.2 0.3 



T/J 

FIG. 2: (a) Calculated J-band for T = 0.07J, 0.23J, and 
0.29J, at a fixed disorder strength a = 0.2J. A one- 
phonon spectral density was used of the form Eq. flU , with 
W (1) = 25 J and a = 3, while the monomer radiative rate was 
set to 70 = 1.5 x 10~ 5 J. (b) Temperature dependence of the 
width A(T) of the calculated J-band (symbols) for three val- 
ues of the disorder strength, a = 0.1J, 0.2J, and 0.3J. Other 
parameters as in (a). The solid, dashed and dotted curves 
represent the corresponding fits, according to Eq. 1291 . 



We next turn to the absorption band calculated ac- 
cording to Eq. I|28[) . neglecting the role of two-phonon 

scattering [W^ = 0]. In Fig. |2fa) this band is plot- 
ted for three temperatures at a fixed disorder strength of 
(7 = 0.2J. For the one-phonon spectral density we used 
the form Eq. (jTTJl with Wq — 25J and a = 3, and took 
7o = 1.5 x 10 _5 J, which is typical for J-aggregates of 
polymethine dyes. Chains of N — 500 molecules were 
considered. The simulated spectra clearly demonstrate 
the thermal broadening, caused by growing homogeneous 
widths of the individual exciton transitions. At low tem- 
perature, the homogeneous broadening is negligible, the 
J-band is inhomogeneous, with a width that is deter- 
mined by the disorder strength. With growing temper- 
ature, the J-band becomes more homogeneous, as is ap- 
parent from the fact that it gets more symmetric. 

In Fig. Hl^b) we plotted by symbols the simulated J- 
band width A(T) as a function of temperature for three 



values of the disorder strength: a = 0.1J, 0.2J, and 0.3J 
[all other parameters were taken as in Fig. Ufa)]. As 
is seen, the A(T) shows a plateau at the value of the 
inhomogeneous width, A(0) = 0.04J, 0.1J, and 0.18J, 
respectively. Beyond these plateaus, A(T) goes up quite 
steeply, reflecting the fact that the homogeneous (dy- 
namic) broadening becomes dominant. To accurately 
extract at low temperatures the small homogeneous con- 
tribution to the total width, we generated up to 4 x 10 s 
disorder realizations. At higher temperatures, this num- 
ber could be restricted to 4000, owing to the reduction 
of the relative fluctuations (cf. Fig. . 

Inspired by the analytically obtained power-laws for 
the one-phonon-assisted dephasing rate as a function of 
temperature for homogeneous aggregates (Sec. lIV'ATl . we 
considered a parametrization of the form 

A(T) = A(0) + alT (1) (T/J) p (29) 

for the total band width in disordered aggregates. It 
turned out that the calculated line widths as a function of 
temperature could be fitted very well by the relation l|29|) 
[curves in Fig. 0T)]. The corresponding fit parameters 
are a = 1.24 and p = 4.16 for a = 0.1J, a = 1.32 and 
p = 4.29 for a = 0.2 J, and a = 1.20 and p = 4.27 
for a = 0.3J. The scaling relation Eq. I)29|) turns out 
to hold over an even wider range of a and Wq values^ 
This implies that, although the J-band is built up from 
a distribution of exciton states with different dephasing 
rates, the total width A(T) may effectively be separated 
in an inhomogeneous width, A(0), and a dynamic con- 
tribution. 

We note that the fitting exponent p is larger than the 
value 3.85 found in the absence of disorder (Sec. IIV Xj . 
This increase results from downward scattering processes 
between optically dominant exiton states, which are pos- 
sible in the presence of disorder, but not for the super- 
radiant state in the homogeneous chain. This claim may 
be substantiated by considering the dephasing rates of 
the lowest exciton state of each disorder realization. We 
numerically generated the average of this quantity, T| , 
from 3 x 10 4 disorder realizations for chains of N = 500 
molecules with a = 0.1J, Wq 1 ' = 25J, and a = 3. This 
average is shown as a function of temperature in Fig. |3| 
(diamonds), together with the dynamic contribution to 
the total J-bandwidth, A(T) - A(0) (solid line), and the 
dephasing rate T^ 1 ' of the superradiant state for a homo- 
geneous chain of the same length (squares) . As is seen, 
the dynamic part A(T) — A(0) has a larger exponent 

p = 4.16 than IY 1 (p — 3.85). It is remarkable, however, 

that r^ 1 ' and display almost identical behavior, at 
least in the relevant region T > A(0), where the homo- 
geneous contribution to the J-bandwidth is noticeable. 
For T <C A(0), where rl 1 ' undergoes significant fluctua- 
tions (see Fig.^l, IV 1 turns out to be much smaller than 
Ii . In the high-temperature regime, these fluctuations 
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are washed out and the two quantities become almost 
identical. 




T/J 

FIG. 3: Temperature dependence of two measures for the ho- 
mogeneous width for disordered aggregates (o = 0.1J) com- 
pared to the dephasing rate of the lowest (superradiant) ex- 
citon state in a homogeneous aggregate (a = 0; squares). For 
the disordered case, we plotted the dynamic part A(T) — A(0) 

of the total J-bandwidth (solid line) and the average value I^ 1 
of the dephasing rate of the lowest exciton state found in each 
disorder realization (diamonds). In all cases we used a chain 
length of N = 500, a monomer radiative rate 70 = 1.5x 10 -5 J, 
and a one-phonon spectral density of the form Eq. lilt with 



W. 



(i) 



from those obtained without fluctuations (i.e., uj = 00: 
triangles). The reason is that at elevated temperatures 
the function aj 3 n(oj) varies slowly on the scale of Q, so 
that the modulating function 1 + sin(27rw/d)) may be re- 
placed by its average value, which equals unity. 




FIG. 4: Temperature dependence of the calculated width 
A(T) using three different models for the one-phonon spec- 
tral density. Diamonds (triangles) correspond to a spectral 
density of the form Eq. 11 U with a = 1 and Wg = 5J 
(a = 3 and W (1) = 25 J), while the squares are the results 



for F m (u)) = W (1) (a;/J) a [l + sin(27ra;/w)] with W U 
and ui — J/6. The various curves are fits to the power-law 
Eq. fZQ. In all cases we used a chain length of TV = 500 and 
a monomer radiative rate 70 = 1.5 x 10~ 5 J. 



25J 



25J and a 



So far, we have only presented numerical results for 
a one-phonon spectral density Eq. (|llfl with the power 
a = 3, which corresponds to a Debye model for the host 
vibrations. To end this subsection, we will address the 
effect of changing the spectral density. First, we consider 
the effect of the value for a. The diamonds in Fig. 0] 
present our results for the temperature dependence of 
the width A(T) obtained for a = 1 and W = 5J. As 
before, we found that these data may be fitted by a sim- 
ple power-law of the form Eq. H29|) (dashed line) . In this 
case we find p = 1.9, which, again, is slightly larger than 
the value 3/2 found from Eq. I|22l) . due to the correc- 
tion of the dispersion relation arising from the long-range 
dipole-dipole interactions and downward scattering pro- 
cesses that contribute to the total J-bandwidth. Clearly, 
these data demonstrate the sensitivity of the temperature 
dependence of the total linewidth to the power a. 

Interestingly, it turns out that while the temperature 
dependence of the J-bandwidth is sensitive to the overall 
frequency scaling of the spectral density, it is not sen- 
sitive to fluctuations of this scaling around an average 
power-law. To demonstrate this, we have considered a 
spectral density of the form T^\uj) = W^Iuj/ J) 3 [l + 
sin(27ru>/(i)], which only on average exhibits an in- 
dependence. The results for the width A(T) obtained 

for Wq 1 ^ — 25 J and uj = J/6 are presented as squares in 
Fig. 2| Remarkably, these results are indistinguishable 



B. Two-phonon-assisted dephasing 

As we have seen in Sec lIII Bl the two-phonon-assisted 

(2) 

dephasing rate TJ, consists of four contributions, one of 
which is elastic (indicated as "pure" , as it is responsible 
for pure dephasing), while the other three are inelastic 
and are indicated as downward (!!), cross (tl), and up- 
ward (It)' In Figs. WL. C )~ Hf); El and|Hlwe present re- 
sults for the statistics of these various contributions for 
disordered aggregates. In all cases, we used chains of 
500 molecules, a disorder strength of a = 0.135J, and 
a two- vibration spectral density of the form Eq. i|ltjfl . 
with a = 3. The statistics are presented for the lowest 
exciton state in each one of 3 x 10 4 randomly generated 
disorder realizations. For this state the downward contri- 
bution vanishes and the two-phonon-assisted dephasing 

rate reads = +Tpure. These three remain- 

ing contributions were calculated using the expressions 
derived in the Appendix. 

From Figs, ^c) andQJd) we observe that at a given 

(2) 

temperature, the relative spread in is much smaller 

(2) 

than that in TW . The reason is that in a two-phonon- 
assisted upward process the exciton in the lowest state 
scatters to more higher-energy states than in a cross pro- 
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(2) 

cess. As a result, fluctuations in are suppressed more 



(2) 

than those in Ti,'. Upon heating, the spread of both 



rl^P and T\V reduces, which has the same explanation as 



.(2) 
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IT 



U 



given for this effect in the case of one-phonon-assisted 
dephasing (Sec. IV A|) . The distribution of r^ re does 
not depend on temperature at all, because, according 
to Eq. I|17fl . the rate r P u re fluctuates exclusively due to 
fluctuations in the inverse participation ratio Xm=i ftn- 
As the latter quantity is subject to large fluctuations, 23 
this also explains the large relative spread in Tp 



0.06 



0.04 



(2) 



0.02 





(2) 

FIG. 6: Temperature dependence of the average values (T^j , 

— (2) — (2) 

r T | , and r pure ) of the three contributions to the two-phonon- 
assisted dephasing rate of the lowest exciton state in disor- 
dered aggregates of N = 500 molecules. Parameters were 
chosen as in Fig. The vertical dotted line shows the cross- 
ing point at T = 0.12J. 



well by a power-law, a(T / J) p . In doing so, we obtained 



= 0.18IU (2) (T/J) 



n8.3 



(30a) 



FIG. 5: As in Fig. Q but now for Fpur e . As is discussed in 
the text, this distribution is solely due to fluctuations in the 
inverse participation ratio of the exciton states and does not 
depend on temperature. 

In Fig. H3 we plotted (on a log- log scale) the tempera- 

—(2) —(2) —(2) 

ture dependence of the mean values I\j , , and r purc , 
obtained from averaging over these rates for the lowest 
exciton states in the simulations discussed above. This 
figure nicely shows the relative importance of the dif- 

_(2) 

ferent dephasing channels. We clearly see that P** <C 

(2) (2) 

T|| ,r pure , i.e., the inelastic channel of dephasing due 
to double phonon absorption is inefficient, at all tem- 
peratures. More importantly, we observe that the cross 
channel of inelastic two-phonon-assited dephasing suc- 
cessfully competes with the pure dephasing contribution. 
At low temperatures, pure dephasing dominates, while at 
higher temperatures the inelastic cross process is more 
important. This correlates well with our findings for 
disorder-free aggregates [see discussion below Eq. l(S7)l]. 
For the disorder strength considered here, the cross-over 
occurs at T Q = 0.12J (w 100 K for J = 600 cm" 1 ). We 
note that when considering two-phonon scattering, one 
often restricts to modelling the elastic process (see, e.g., 
Ref . WO) . From the above we see that this is not justified 
at elevated temperatures. 

Despite the fact that not all three curves in Fig. are 
exactly straight lines (deviations occur in the unimpor- 
tant low-temperature part), they all can be fitted very 



r}f = 9.01IU (2) (T/J) 7 ' 9 , (30b) 



rjfl = 1.88IU (2) (T/J) 7 . (30c) 

We recall that the exponent p = 7 in the last formula is 
an exact result for a = 3 and lo c 3> T, as was argued at 
the end of Sec. lIV El already. As in the case of one-phonon 
scattering, we see that the inelastic two-phonon dephas- 
ing rates exhibit a steeper temperature dependence than 
for the homogeneous aggregate with nearest-neighbor in- 
teractions [Eq. 

VI. COMPARISON TO EXPERIMENT AND 
DISCUSSION 

It is of interest to see to what extent the model we 
presented here is able to explain the temperature depen- 
dence of the homogeneous broadening in molecular ag- 
gregates. As mentioned in the Introduction, Renge and 
Wild 11 found that the total J-bandwidth A(T) of PIC-C1 
and PIC-F over a wide temperature range (from 10 K to 
300 K) follows a power-law scaling as in Eq. 129[) . Al- 
though the power reported by these authors (p — 3.4) is 
smaller than the ones we derived in Sec. IV Al from direct 
comparison to the experimental data we have found that 
a model of one-phonon scattering with a spectral density 
given by Ea. l|ll|l with a — 3 and oj c — > oo, yields an ex- 
cellent quantitative explanation of the experiments over 
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the entire temperature range, both for the shape and the 
width of the J-band. The same turns out to be true for 
the J-bandwidth of PIC-Br, measured between 1.5 K and 
180 K>£ In all these fits, a and Wo were the only two free 
parameters that could be adjusted to optimize compari- 
son to experiment. Details will be published elsewhere^ 
together with a fit of the much debatecS^^^ 2 ^ tem- 
perature dependence of the fluorescence lifetime of these 
aggregates. 

Here, we present an explicit comparison to the hole- 
burning data reported by Hirschmann and Friedrich^S 
Using this technique, they measured the homogeneous 
width of the exciton states in the center of the J-band 
for PIC-I over the temperature range 350 mK to 80 K. 
Their data for the holewidth T are reproduced as trian- 
gles in Fig. The solid line shows our fit to these data, 
obtained by simulating disordered chains of N — 250 
molecules with a one-phonon spectral density of the form 
Eq. Hll|l with a — 3 and lu c — > oo. The resonant in- 
teraction strength and the monomer radiative rate were 
chosen at the accepted values of J = 600 cm -1 and 
7o = 1.5 x 10~ 5 J = 2.7 x 10 8 s _1 , respectively. Thus, 
the only free parameters were the disorder strength a and 
scattering strength W^p . First, we fixed the value of a by 
fitting the low-temperature (4 K) absorption spectrum, 
where the homogeneous broadening may be neglected. 
This yielded a = 0.21J. Next, W (1) was adjusted such 
that the measured growth of the hole width was repro- 
duced in an optimal way. Thus, we found Wq = 180J. 
For each given temperature the simulated hole width was 
obtained as the average of the dephasing rate of the exci- 
tons found in an interval of 0.06 cm -1 around the center 
of the simulated J-band at that temperature. 
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FIG. 7: Hole width T as a function of temperature measured 
for aggregates of PIC-I after pumping in the center of the 
J-band (triangles) i— The solid line is our fit, assuming one- 
quantum scattering on acoustic phonons in the host. Model 
parameters are discussed in the text. 

Taking into account the fair amount of scatter in the 
experimental data, we conclude from Fig. {7\ that our 
model yields a good fit to the measurements. As men- 



tioned above, the same holds for the J-bandwidth in ag- 
gregates of PIC-C1, PIC-F, and PIC-Br. This yields valu- 
able information about the dominant mechanism of de- 
phasing in these materials. We conclude that this mech- 
anism is one-phonon-assisted scattering of excitons on 
vibrations in the host characterized by a spectral density 
which (on average) scales as u 3 ; this scaling yields a nat- 
ural explanation of the power-law thermal broadening of 
the J-band found in various experiments. We stress that 
it is impossible to fit the experimental data with a spec- 
tral density that is constant or scales linearly with u>, 
as that yields considerably different power laws for the 
width [cf. Fig. 0|- The ui 3 scaling of the spectral den- 
sity needed to fit the experiments strongly suggests that 
acoustic phonons dominate the scattering process. Thus, 
the spectral width measured over a broad temperature 
range is an excellent probe for the scattering mechanism. 

It is appropriate to comment on the value of Wq 1 ^ ob- 
tained from our fit, which seems to be very large. It 
should be kept in mind that Wq 1 ^ is a phenomenological 
scattering strength, which combines several microscopic 
material properties [see discussion below Eq. fljlXp]. Most 
importantly, the value found here is consistent with a per- 
turbative treatment of the scattering process: the scat- 
tering rates between the optically dominant states ob- 
tained from it turned out to be much smaller than their 
energy separation. 

We finally address an alternative mechanism of de- 
phasing, namely scattering on local vibrations belong- 
ing to the aggregate. This has been suggested by sev- 
eral authors based on activation-law fits of the mea- 
sured homogeneous contribution to the J-bandwidthii£ 
To get an estimate whether this is a reasonable mecha- 
nism, let us neglect the disorder and make the nearest- 
neighbor approximation for the resonance interactions 
J nm . Then, E v and tp url are given by Eq. I|19|) . Close 
to the lower exciton band edge, the region of our inter- 
est, E u = -2J + Jk 2 v 2 /{N + l) 2 . Furthermore, let us 
parameterize the spectral function of a local vibration of 
frequency ujq as T(u>) = 2t:Vq5{uj — u)q), where Vq is the 
coupling constant to the excitons. We are interested in 
the dephasing rate Ti of the superradiant state \v = 1). 
Using the above simpifications and replacing in Eq. I|12|) 
the summation over exciton states by an integration, one 
easily arrives at 



ri 



2 n(cj ) 
(Juo) 1 / 2 ' 



(31) 



Comparing this result to the activation law 
6exp(— ljq/T), used in Refs. and to fit the 
experimental data, one obtains U 2 = b(Jujo) 1 / 2 . Sub- 
stituting the values b = 3000 cm -1 and loq = 330 cm -1 
from Ref. [ToL and J = 600 cm -1 , we obtain as estimate 
for the exciton-phonon coupling Vq » 1100 cm -1 . This 
is an enormously large number. In particular, the Stokes 
losses S = Vq/ujq ~ 3700 cm -1 turn out to be much 
larger than the resonant interaction J = 600 cm -1 . 
Under these conditions, a strong exciton self-trapping is 
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to be expected, resulting in a reduction of the exciton 
bandwidth by a factor of exp(— S/w ) "C 1^ Coherent 
motion of the exciton, even over a few molecules, is 
then hardly possible: any disorder will destroy it. 
This observation is not consistent with the widely 
accepted excitonic nature of the aggregate excited 
states, as corroborated by many optical and transport 



measurements 



3.45 



VII. CONCLUDING REMARKS 

In this paper we presented a theoretical study of the 
temperature dependence of the exciton dephasing rate in 
linear J-aggregates and the resulting width of the total 
absorption band (the J-band). As dephasing mechanism 
we considered scattering of the excitons on vibrations of 
the host matrix, taking into account both one- and two- 
vibration scattering. The excitons were obtained from 
numerical diagonalization of a Frenkel exciton Hamilto- 
nian with energy disorder and their dephasing rates were 
subsequently calculated using a perturbative treatment 
of the exciton- vibration interaction (Fermi Golden Rule) . 

In the absence of disorder, the lowest (superradiant) 
exciton state dominates the absorption spectrum. As a 
result, the dephasing rate of this state directly gives the 
homogeneous width of the J-band. We analytically cal- 
culated the temperature dependence of this homogeneous 
width for both one- and two- vibration scattering, assum- 
ing a Debye-like model for the host vibration density of 
states. It turned out that in all cases the homogeneous 
width obeys a power-law as a function of temperature, 
with the value of the exponent depending on the shape of 
the low-energy part of the vibronic spectrum CSec. irVl . 

In the presence of disorder the optically dominant ex- 
citon states still reside close to the bottom of the band, 
but their energies are now spread and their wave func- 
tions become localized on finite segments of the chain. 
We have found that the various one- and two-phonon- 
induced contributions to the dephasing rate undergo sig- 
nificant fluctuations, because the exciton energies and 
overlap integrals vary considerably from one disorder re- 
alization to the other (Sec. El . As a consequence, one 
cannot use the dephasing rates of individual states to 
characterize the homogeneous broadening of the J-band. 
Instead, we simulated the total J-band and showed that 
its width effectively separates in an inhomogeneous (zero- 
temperature) contribution and a dynamic (homogeneous) 
part. The latter scales with temperature according to a 
power law [Eq. (|29|l ] , with an exponent that is somewhat 
larger than the one found for the homogeneous broad- 
ening in the absence of disorder. We also showed that 
amongst the two-vibration scattering processes, inelas- 
tic channels will at elevated temperatures dominate the 
usually considered pure-dephasing contribution. 

Finally, from comparison to absorption and hole- 
burning experiments fSec. IVljl . we found that the domi- 
nant mechanism of dephasing for J-aggregates lies in one- 



phonon scattering of excitons on vibrations of the host 
matrix, characterized by a spectral density which (on av- 
erage) scales like the third power of the phonon frequency. 
This suggests that acoustic phonons of the host play an 
important role in the scattering process. All tempera- 
ture dependent data available to date, are consistent with 
this picture. By contrast, we argued that the previously 
suggested mechanism of scattering on local vibrations of 
the aggregate, leading to an activated thermal behavior, 
is not consistent with the overwhelming amount of ev- 
idence that the optical excitations in J-aggregates have 
an excitonic character. 



APPENDIX A: TWO-PHONON SCATTERING 
RATES 



In this Appendix we present expressions for the 

(2) 

two-phonon-assisted scattering rate W^J , starting from 
Eq. I|15|) . After substituting the two- vibration spectral 
density given in Eq. (|16f) and performing several alge- 
braic manipulations, we arrive at 



2^ VpnVvn 



(2) 

where the temperature dependent functions (oj^), 

(2) (2) 

F^y' {ujfiv), and Fi^u^v) distinguish between the scat- 
tering processes in which two phonons are emitted (||), 
one is absorbed and one is emitted (tl), and two phonons 
are absorbed (TT)- They are given by 



(Al) 



JT ,u; c /T 

dx I dy 



fll(x,y)8(^- + x + y) , (A2a) 



u} a /T ruJ c /T 

dx i 

o JO 



F^f^^) = I dx I dy 



fn(^y)s(^f-x + y) , (A2b) 



•Ft? (<*V) = / dx j dy 
Jo Jo 



x f^(x,y)d(^f-x-y) , (A2c) 

where, after changing to dimensionless integration vari- 
ables x = ojq/T and y = u q i /T, we introduced the auxil- 
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iary functions 



fll{x,y) = 

fn(x,y) = 
frtfav) = 



x p e x 


y p e y 


e x _ 


1 


e» - 1 


x p 




y p e y 


e x — 


1 


ev - 1 


x p 




yP 


e x _ 


1 


e» - 1 



(A3a) 
(A3b) 
(A3c) 



We note that for /i = v = 0), the only non vanishing 
term is F^(0), which describes the elastic channel of 
scattering (pure dephasing). 

We now further analyze the three contributions to the 

(2) 

scattering rate. As T^{u>^ v ) describes the emission of 
two vibrational quanta, we have 10^ < 0. Thus, per- 
forming the y integration, we obtain 

F[f(ov)= / ^ dxf u (x,-^-x) ,(A4a) 



T 



if — < uj c , and 



*lf(<*V) = / dxf u (x,-^-- x) 



if lo c < —tjj^u < 2oj c , while (u)^) = otherwise. 

Next, the contribution which involves the emission and 
absorption of one vibrational quantum, is given by 



f T ( fK„)= dxf n (x,-^+x) , (A5a) 



if < Ul^ < LO c . If < -ui^ < u> c , 
n(w c +w fl „)/T 



J o dxf n (x,-^+x) (A5b) 



and F^{uj^ v ) = otherwise. 

Finally, the contribution that results from the absorp- 
tion of two vibrational quanta (w M „ > 0), yields 



r ^ x 

T 



) , (A6a) 



if < uOn U < uj c , and 



7 (2) 

TT 



- dx/ TT (x, ^ - *) , (A6b) 



( 2") 

(A4b) if o> c < < 2cj c , while F| T ; (cj My ) = otherwise. 



T1 



* On leave from S.I. Vavilov State Optical Institute, Birzhe- 
vaya Liniya 12, 199034 Saint-Petersburg, Russia. 
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